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Abstract 

We compute the full x - dependence of the proton and neutron 
spin structure functions in the MIT bag model, including the effect of 
gluon exchange and the meson cloud. Impressive agreement is found 
for x larger than 0.1, where polarised gluons are not expected to play 
a significant role. 
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At the present moment there is an impressive collection of high energy 
data on nucleon structure functions that defy a complete understanding. To 
quote a few, we mention the NMC |TJ measurement of the Gottfried sum rule, 
the NA51 p[ measurement of the sea asymmetry and the (already classical) 
problem with the polarized structure functions as measured by EMC/SLAC 
[3 - 6]. In the particular case of the polarized structure functions, it now 
seems that one has reached the end of the "spin crisis" and the beginning of 
the spin problem @, 0|. As it is unlikely that perturbative QCD corrections 
alone are enough to solve the problem, we should examine the possibility 
that at least some of the discrepancy arises in the non-perturbative regime. 
A common approach has been to calculate the desired quantities at some low 
scale and then evolve them, using the QCD renormalization group equations, 
to the experimental scale. Unless the model has been specifically derived from 
QCD, the starting scale at which it best approximates the structure of the 
nucleon is unknown. It is usually adjusted to fit the experimental parton 
distributions as well as possible. Of course, one might expect that some 
models would give a better description of the data than others and there is 
considerable interest in using DIS to help choose amongst models |§. 

Although there are a variety of possibilities to calculate quark distibutions 
at the starting scale, like non-relativistic quark models |§ or vertex functions 
1(J , here we shall work with the MIT bag. This choice is based on its success 



in low energy physics, its simplicity and the insights it has already yielded 
in connecting low and high energy data fll], [L2| . Of course we do not expect 
that the bag model alone contains all the physics we need because we know 
that mesons should play a role as well. For instance, we know that chiral 
symmetry in the bag is restored through a meson cloud. Moreover, the meson 
cloud is already known to give important corrections to some sum rules [14 
- 19] 

As important as it is to correct the integrated quark distributions, the 
study of the effect of the meson cloud on the x dependence of the quark 
distributions themselves is at least as important [21 - 24]. In this letter we 
calculate the x dependence of the polarized and unpolarized quark distribu- 
tions in the proton in a bag model dressed by mesons. We perform the cal- 
culations in next-to-leading order (NLO) QCD. Moreover, we present for the 
first time the x dependence of the combined effect of the N - A interference 



term |2T| and one gluon exchange | 22| . The interference term is important 



because, as noticed before |I4]|, in the Ellis- Jaffe sum rule it cancels part of 
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the reduction coming from pions. As we shall see, the interference terms 
have a determinant role in the shape of the polarized quark distribution of 
the neutron. 

We introduce mesons in the model through the Sullivan process [p5|1 , 
simply noting that there are unresolved questions about the model, partic- 
ularly the validity of the impulse approximation p6 |. The modern study 
of the mesonic contribution via convolution started with Thomas [27] and 
was later extended to the study of structure functions by the Adelaide group 
| fiG| . |24 , 28] . A problem remaining in these calculations is the freedom for the 



value of the cut-off in the form factor. To avoid this problem, we shall follow 



the approach of the Julich group |L9|, P9[ , where the cut-off is fixed through 
the use of high energy pp data. 

The basic hypothesis in this sort of model is that the physical nucleon 
wave function (in the infinite momentum frame) can be written as a super- 
position of a few Fock states: 



\N) bare + Y^ [ dx f d 2 k ± cj) BM (x, k±)\B(x, k±),M(l - x, -k±)) 

BM J ° J 



\N) phys = Z 1 ' 2 

The wave function renormalization factor 

1 + [ dx [ d 2 £a|(W2,A; ± )| s 

BM J ° J 



(2) 



measures the probability that the physical nucleon contains a bare nucleon. 
The Fock states used in our calculation involve the low mass particles which 
are important to describe nucleon properties, namely the nucleon SU (3) octet 
(N, A, E) and decuplet (A, £*) and the first pseudoscalar and vector meson 
octet (tt, K, p, u, K*). We have included the hyperon-kaon contributions 
for completeness, but their actual size is very small (~ 2%) Jl9|, 20, 28 1. 



The Fock state expansion Eq.([l|) has consequences for the structure func- 
tion of the nucleon. Due to the presence of baryon-meson Fock states, the 
virtual photon can scatter either on the nucleon core or on the meson-baryon 
system. Formally, the quark distribution q(x) of the nucleon is given by: 



q(x) 



N 



[X 



+ J2 ($ M q(x) + s B g(x 



BM 



(3) 
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The contributions from the virtual meson and baryon can be written as the 
convolution of the meson (baryon) structure function with its longitudinal 
momentum distribution in the nucleon: 



S M q(x) 



fMB(y)q M ( ~) ^ and 5 B q(x) = £ f BM (y)q B ^ dy 



m) y 

where f MB and f BM are given by 



(4) 



Jmb(x) = J d 2 k ± \(j) BM (l-x, -k ± )\ 2 , and f B M{.x)=Jd 2 k 1 _\(f) BM (x,k A _)\ 2 . 

(5) 

In order to conserve charge and momentum, we have the following relation: 

/bm(1 - v) = fMB{y)- (6) 
The functions f B M(y) and fMB^y) can be calculated using time ordered per- 



turbation theory in the infinite momentum frame [2g, ^9], pi 
forms for fsuiv) and /ms(z/); can be found in Refs. 



The analytic 



mm. 

In practical calculations we need more information, namely the various 
coupling constants and the vertex form factors G B m- The coupling constants 
can be extracted from scattering experiments and are rather well known [^] . 
For the vertex form factor we use an exponential parametrization: 



G 



BM 



exp 



1 



9A 2 



{m 2 N -M 2 BM (y,k 2 ± )) 



(7) 



with M 2 MB 



(y,k 2 



m 2 R +k 2 . 



+ 



i . The cut-off parameters Abm can be 
estimated, using one boson exchange models, from n, A and A ++ production 
in high energy pp scattering (for details see Ref. Hl9|). They were found to be 
Atvtt = A/vTp = 1.08 GeV and Aatt = Aa p = 0.98 GeV. The procedure used 
to obtain the cut-off parameters may be questioned because the assumption 
of single meson exchange being responsible for the process is fairly simple. 
However, it certainly gives an upper bound for the values of the various A s. 

At this point we are still left without the input distributions in equa- 
tions (|3p and (H). For the quark distribution in the pion we use a recent 
parameterization by Sutton et al. 133]. By using SU(3) symmetry, the quark 



4 



distributions of all the mesons can then be obtained. For the bare quark 
distribution of the baryons we use the bag model calculation of the Ade- 
laide group [JTTJ, |53[. The main advantage of this method is that it ensures 
energy-momentum conservation and hence the correct support of the quark 
distributions. It has also been successful in describing data at high momen- 
tum transfer, Q 2 , both in leading order [11] and in next-to-leading order 
[|T^] . Mesonic corrections lead to further improvement in the predictions of 
the model and will, among other things, reduce the amount of evolution 
needed [341. The form of the quark distribution is given by |! I 



in 



X \Pn\d\p n \^^\^iPn)\ 2 . (8) 

i[M2(l-n)2-M2]/2M(l- I ) 103(0)1 

Here is the spin- flavor part of the wave function of the initial state (at 
rest), Pf, m makes the projection onto flavor / and spin projection m, M n 
is the mass of the intermediate state and ip the Fourier transform of the 
quark wavefunction. Equation (§) gives the two quark contribution to the 
total quark distribution which dominates at intermediate and large x. We 
should also account for the contributions coming from four quarks in the 
intermediate state, but for simplicity we shall mimic this contribution by a 
term of the form (1 — x) 7 , properly normalized (further discussion on this 



subject can be found in Ref. |TT 



We use the MRS parametrization [35] of the unpolarized structure func- 



tions to fix the parameters of the model; e.g., the radius of the bag, the 
average mass of the spin scalar (M s ) and spin vector (M v ) diquark in the 
intermediate state and the low scale, fi 2 , at which the model is supposed to 
be valid. In Fig. [I] we show the total valence distribution computed in NL0& 
in the MS scheme for the bag dressed with mesons. Very good agreement 
with the MRS parametrization is found for /i 2 = 0.165 GeV 2 , R = 0.8 fm, 
M s = 0.65 GeV and M v = 0.85 GeV. For comparison, we also show the bare 
bag calculated with the same set of parameters. The bonus of the NLO fit 
is that it provides a decrease of a s at \j? from ~ 1.52 in a LO fit to ~ 0.6 in 
NLO. This is also a drop in comparison with the value a s ~ 0.77 found in 



Details for NLO calculation in the bag can be found in 
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a NLO fit without mesons |12[]. For completeness, we note that these values 
of the coupling constant were found using Aqcd = 0.2 GeV and three active 
flavors. 

After we have fixed all parameters in the unpolarized deep inelastic scat- 
tering sector, we can explore the consequences for the polarized sector. The 
calculation of effects due to the presence of higher BM Fock states for the 
gi(x) structure function is similar to those for q(x) given in Eq.(^). The 
contribution from the scattering on the recoil baryon is 



-l 



' x\ dy 



5 B Aq(x) = J2[ d BM (y)Aq B [ - ] — , 

m J * \yj y 



(9) 



where ds{y) is the polarized, longitudinal momentum distribution. It can be 
calculated using the same techniques as for the unpolarized case (for details 
see HT9| ). The mesonic contribution vanishes because of the pseudoscalar 



character of the pion. The main difference from the unpolarized structure 
function is the presence of the N- A interference term [[H]] , which can also be 



(10) 



written as a convolution ||29|| : 

5 mt Aq(x)= CdUv)^ (-)-■ 
J* \yj y 

The necessary polarized splitting functions dsu and d int can be found in the 
Appendix. 

Combining all these contributions, gi for the nucleon is given by: 



gPhys 



gbare _|_ V' I \Li ,.B 



BM 



(5 B g* + 5g\ nt 



(11) 



We now apply Eqs. (|) - (JTTj) to calculate g\ (x). In the meson sector 
only the pseudoscalars are included. The vector mesons are omitted because 
we are not aware of any model to extract their polarized quark distribution 
functions. Once again we need a model for the bare quark distributions in 
the nucleon and, as before, we use the MIT bag model. For the polarized 
case we need to specify the spin-flavor part of the wave function. We shall 
use the usual SU (6) wave function in which case, for the interference terms, 
we have: 



(p T ,nV|A + ' T ,A°' T ) = (^n^lA+'^A ^) 



V2 



(12) 
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(p T ,nV|A +,T ,A°* T > = (pWl^lAH A°' T ) = " 




We notice that in the matrix elements fljj|) and (|l~3|) only mixed symmetric 
terms contribute and, as a consequence, the intermediate state always forms 
a spin vector. 

The first moment of the polarized structure function for the proton, g\{x), 
is expressed in NLO as: 



iWMi+iK-^BO-^ (14) 

with g a and g% nonsinglet distributions and go a singlet distribution. How- 
ever, the full singlet anomalous dimensions (for any moment) for polarized 
scattering in NLO are still not known and, because of that, it is not possible 
to calculate the x dependence of go in NLO. Faced with this problem, we 
decided to take the following two approaches to the evolution of g^: 

• In the first approach, case (a), we evolve g a and g$ as nonsinglet in NLO 
and evolve go as a singlet in LO so that go does not pick up the (1 — - 1 ) 
correction and the whole structure function is overestimated. 

• In the second approach, case (b), we treat g a , g$ and go as nonsinglet 
combinations and evolve them in NLO. In this case, go picks up a correction 
of the form (1 — - ) such that the corrections due to the NLO evolution 
are overstimated and the structure function is underestimated. The actual 
curve must be somewhere between the two approaches. 

In Fig. ^ we show the EMC and earlier SLAC data for xgf(x) together 
with the bare bag, the bag plus mesons but without the N-A mixing terms 
and the bag plus mesons plus mixing terms. We stress that the parameters 
for the bare bag differ significantly from those used in earlier calculations 



II]. The meson cloud lowers the bag model prediction over the entire range 



of x, in accordance with earlier estimates for the Ellis- Jaffe sum rule in the 



bag []14] . This is because some of the spin of the nucleon is carried as angular 
momentum by the mesons. The actual value of the calculated Ellis- Jaffe 
sum rule at 10 Gey 2 drops from 0.209 in the bare bag to ~ 0.173 in the bag 
plus mesons for case (a) and to ~ 0.169 for case (b). The actual value of 
the sum rule in the bag model plus mesons, calculated using Eq. ( |H| ) as it 
stands, is ~ 0.171, supporting our claim that the full NLO prediction for the 
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x dependence of gi p (x) in the present model is somewhere between the results 
for case (a) and (b) of Fig. ||. The inclusion of vector mesons is expected to 
reduce the value of the sum rule somewhat more. 

The x dependence for the model, compared with the SLAC E143 data || 
for <7i(x), is shown in Fig. [3], and is quite impressive. This data was taken 
at an average 3 GeV 2 and has smaller error bars than earlier experiments. 
Comparison with this set is also a good test of our model once we have to 
move to a different Q 2 . The resulting agreement between theoretical and 
experimental values for g\(x) is inspiring and provides some confidence in 
the model. At this value of Q 2 the calculated value for the Ellis- Jaffe sum 
rule is ~ 0.171 for case(a) and ~ 0.166 for case (b). 

These results shed light on how the spin in the physical proton is shared, 
suggesting strongly that mesons are responsible for part of the dilution of the 
spin. In general, we can say that the agreement between the data and the 
theoretical calculation is very impressive and that further corrections might 
well bring the entire curve within experimental errors. For instance, we know 
that in the axial gauge the axial anomaly is related to the polarized gluon 
distribution |36 |. In this context, we think that a reasonable polarized gluon 
distribution [37] could bring the curve down in the region x < 0.4. 

We would also like to call attention to the role of the N-A mixing term. 
As said before, it tends to increase the Ellis- Jaffe sum rule. In fact, if the 
mixing terms were absent, the new value of the sum rule at Q 2 = 10 GeV 2 
would be ~ 0.159 for case (a) and ~ 0.155 for case (b). For the Q 2 compatible 
with the SLAC E143 data, these values would be reduced to ~ 0.158 for case 
(a) and to ~ 0.152 for case(b). Figures 2 and 3 tell us that the rise of g\{x) 
due to the mixing terms is confined to the region x < 0.3. This is because 
for the mixing term, the mass of the intermediate state is always M v , and 
then the contribution is isolated at smaller x when compared with the other 
contributions (similar to the down quark distribution in the bag 

The most interesting effect associated with the mixing terms can be ob- 
served in Fig. |], where is shown. Although in gi(x) the effect of the 
mixing terms in the x distribution is not too dramatic because one is adding 
a small number to a large number, in g™(x) the effect is relatively large be- 
cause one is subtracting from numbers near to zero. In fact, for g™(x) the 
mixing terms are essential to give to the theoretical curve the shape of the 
experimental data as measured by the SLAC E142 || experiment. The first 
moment of the calculated polarized distribution of the neutron turns out to 
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be ~ 0.004 for case (a) and ~ —0.003 for case (b). It is also worth noticing 
that there is a consistency between the calculation of gi p (x) and gi n {x) - the 
same calculation that fits the unpolarized data also makes a good prediction 
for both gip(x) and g\ n (x). 

We would like to thank S. Bass, W. Melnitchouk, G. Piller, A. W. 
Schreiber and J. Speth for valuable discussions during the development of 
this work. One of the authors (H.H.) would like to thank the German Aca- 
demic Exchange Service (DAAD) for financial support for his visit to Ade- 
laide. This work was supported by the Australian Research Council, by the 
DAAD and by CAPES (Brazil). 



A Appendix 



y 2 (l-y) [m%-M^(y,kl)] 



2 ' 



(A.l) 




x 



[(ym N + m A ) 2 + k 2 ^[(y 2 m 2 N - m 2 A ) 2 + &ym N m A k\ 
W-ML(y,fci)]2 



^A.2) 



and 




x 



[ -m N (l - y)(ym N + m A ) 2 (ym N - m A ) 
I [™%-MUy,k 2 ± )][m 2 N -M 2 N7T (y,k 2 ± )} 



(2m A + (3y — 2)rriNm A — ym 2 N )k 2 ± — k\ 



(A.3) 



[m 2 N -MUy,kl)][m 2 N -M 2 N ^kl)\ 
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Figure 1: Total valence distribution in the bag and in the bag dressed with 



mesons compared with the MRS [3l| parametrization of the data in the MS 
scheme. The quark distributions are evolved in nest-to-leading-order QCD. 
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Figure 2: Polarized quark distribution of the proton as measured by the EMC 
collaboration || against theoretical predictions for a bare bag, a bag with 
mesons without mixing terms and a bag with mesons and mixing terms. 
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Figure 3: Polarized quark distribution of the proton as measured by the 
SLAC-E143 experiment against theoretical predictions for a bare bag, a 
bag with mesons without mixing terms and a bag with mesons and mixing 
terms. 
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Figure 4: Polarized quark distribution of the neutron as measured by the 
SLAC-E142 experiment || against theoretical predictions for a bare bag, a 
bag with mesons without mixing terms and a bag with mesons and mixing 
terms. 
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